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Abstract 

Instanton - antiinstanton pair is considered as a source of singularity at the Borel 
plane for the ground state energy of anharmonic oscillator. The problem of defin- 
ing the short range instanton - antiinstanton interaction reduces to calculation of a 
smooth part of the Borel function, which cannot be found without explicit calculation 
of several terms of ordinary perturbation theory. On the other hand, the large order 
terms of perturbative expansion are dominated by large fluctuations in the functional 
integral like well separated instanton and antiinstanton. 

The preasymptotics (~ 1/n) of large order perturbation theory contribution to 
the ground state energy of anharmonic oscillator was found analytically. To this end 
the subleading long range asymptotics of the classical instanton - antiinstanton inter- 
action, the one - loop quantum contribution to instanton - antiinstanton interaction 
and the second quantum correction to a single instanton density were considered. 
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1 Introduction 



In this note we would like to consider how the large fluctuations in the Euclidean functional 
integral contribute to the ground state energy of double well anharmonic oscillator 

Our practical aim will be the analytical derivation of the ~ 1/n correction to the instanton 
- antiinstanton pair induced asymptotics of the perturbative expansion: 

E = T,K9 2n ,E n = - 6 —^(l-^- + ...)• (2) 

After the works of L.N.Lipatov [0 it is generally recognized, that classical paths in 
the Feynman functional integral contribute significantly to the asymptotics of the ordinary 
perturbation theory. However the attraction between instanton and antiinstanton makes 
the high order estimates for the double well oscillator (fj) somewhat more complicated. 
The usual trick to avoid this problem is the analytical continuation of the instanton - 
antiinstanton contribution over the coupling constant ||. 

In present paper we prefer not to make the ambiguous analytical continuation at the 
intermediate steps of approximate calculation. Following t'Hooft || (see also JU) we trans- 
form the functional integral to the Borel type integral by considering the action as a col- 
lective variable. In this approach the instanton - antiinstanton pair manifests itself as the 
singularity of the Borel transform. The problem of short range instanton - antiinstanton 
interaction reduces to that of calculation the shape of Borel transform far from singular 
point. On the other hand the smooth part of the Borel function is sensitive mostly to a 
first few terms of perturbative expansion. Thus not only instanton helps one to find the 
asymptotics of perturbation theory, but also the precise calculation of several terms of per- 
turbative expansion allows to avoid the problem of short range instanton - antiinstanton 
attraction. 

The type of singularity at the Borel plane is determined by interaction between the 
instanton and antiinstanton. Therefore one have to go beyond the dilute gas approximation 
in order to find the asymptotics of perturbative expansion. Up to now all authors have 
taken into account only the leading ~ e~ T / g 2 term of the classical instanton - antiinstanton 
interaction (T is the distance between the pseudoparticles). In order to find the ~ 1/n 
correction to @ we have to consider the subleading ~ Te~ 2T /g 2 correction to classical 
instanton - antiinstanton interaction, the one loop quantum correction ~ Te~ T to the 
interaction and the second quantum correction ~ g 2 to the single instanton density. 
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2 Asymptotics of the perturbative expansion and non- 
perturbative effects 



The ground state energy for the Hamiltonian (|IJ) is given by the functional integral 
E(g) = lim — - In Tre~ HL , Tre~ HL = N f Dx{t) e~ s , 



S = j\dt(±x 2 + \x\l-gxf) . (3) 
Instanton - the well known classical path reads: 

g{e + 1) vg 

This path connects two distinct vacua of the action (]3|) x = and x = 1/g. We want 
to find the configurations, which make a large contribution to the high order terms of 
perturbation theory and therefore may be considered as a fluctuation (although large) 
around the trivial vacuum x = 0. The example of such topologically trivial fluctuation is 
the instanton - antiinstanton pair. The pair is not the exact solution of the equation of 
motion and therefore we have to understand, what the specific configuration may be called 
the instanton - antiinstanton pair? But to a first approximation the interaction between 
pseudoparticles is model independent for a wide class of reasonable configurations (see e.g. 

0) 
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S[x I - A } = -(l/Z-2e- T + •••) . (5) 



,2 
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The instanton - antiinstanton contribution to Tre~ HL is given by: 

(Tre-HLy-A „ L j ^L exp{ _ J_ (1/3 _ 2e ~ T + •••)}• (6) 
J ^9 9 

Factor L appears after integration over the center of the pair and (7rg 2 ) -1 comes from the 
square of single instanton density (see 0). In order to find the ground state energy one 
have to sum over the "dilute gas of instanton - antiinstanton pairs". Hence the energy 
reads: 

E(g) = -( — 2 exp{-- 2 (l/3 - 2e~ T )}. (7) 
J ng 2 g 2 

This integral diverges at large T. The divergency will disappear if for large T one replaces 
the summation over dilute gas of pairs by summation over dilute gas of individual pseu- 
doparticles. Nevertheless even the expression (0) carries quite sufficient information about 
high orders of perturbation theory. Consider new variable x = g 2 S, where S — is a classical 
action of the instanton - antiinstanton pair. The value of x varies from x = (the classical 
vacuum) to x = 1/3 and for widely separated pseudoparticles x = 1/3 — 2e~ T + ... . Now 
formula ([I]) takes the form: 

E= f 1/3 F(x)e-^ 2 — 2 , F(x)*- 1 (i + (gg)) + ... . (8) 
Jo Tig 2 [1/6 — x) 
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About the function F(x) we know only that it has a pole at x = 1/3. However just the 
singularities of F(x) at complex x plane determine the behaviour of series F(x) = F n x n at 
large n and thus the asymptotics of perturbative expansion of E(g). The only assumption 
that pole at x = 1/3 (|8|) is the closest to origin singularity of F(x) allows one to find 
the well known |7j leading term of the asymptotic expansion (H). Below we will find the 
weaker (logarithmic) singularity of function F(x) at the same point x = 1/3. There is no 
direct way to calculate the smooth function F(x) at x < 1/3, but after one have calculated 
precisely several terms of E(g) perturbative expansion, the smooth part of F(x) carries 
now new information. 

At very large T the dilute gas of many pseudoparticles is to be considered. In or- 
der to eliminate the double counting one have to subtract from (|7]) the contribution of 
noninteracting instanton and antiinstanton: 

/ r T=0 ° dT sm f T =°° dT e- 1/6g2 

E ± = -( —e- S(T) - / — e" 1/39 =F ^ ■ (9) 

\Js(T)=o irg z Jt=o irg z / g^Jix 

Here the last term describes the nonperturbative splitting of the ground state due to tun- 
neling through the barrier. Two signs correspond to the states having different parity with 
respect to transformation x — > 1/g — x. After passing to variable x = g 2 S the eq. © takes 
the form: 

/ fi/3 e - x '9 2 - e- 1 / 3 ^ dx e- 1 / 3;?2 ln(6)\ e- 1 / 6 ^ , . 

E± = ~\ — m 2 2 — J t — t^- ( 10 ) 

\Jo 1/3 — x ng A ng z J g^Jix 

This expression correctly accounts for the high order terms of perturbation theory and the 
main nonperturbative effects caused by single instanton or instanton - antiinstanton pair. 
However the first terms of perturbative expansion are completely wrong in ([TIT) . In order 
to make ([II]) quantitatively reliable one have to calculate precisely the iV first terms of the 
perturbation theory and subtract from ( [TOD the N first terms of its expansion in series in 

9 2 : 



J2n 



" f i-i/3 e - x 'v - e-W dx e- 1 /3ff 2 z n (6) 3 n+1 (n)! , 

E ± = E n - < — 2^ 9 

1/3 -x irg 2 ng 2 ^ u 

p-l/6 9 2 

T—y^il + E^K) ■ (11) 



g^/n 



n=l 



Here the last term ~ exp(— l/6g 2 ) formally is the largest nonperturbative correction, but 
this term cancels if one considers the "center of band" E cb = \{E + + EJ). The terms of 
series (0) begin to grow up at n > l/3g 2 . Respectively the approximate expression (|TT|) 
reaches the best accuracy at iV = l/3g 2 . Calculating the new terms of perturbation theory 
at iV > l/3g 2 one should only decrease the accuracy of total results. 



In addition to the "center of band" the formula (JlTT) allows to calculate also the indi- 
vidual energy levels. One have to calculate only a lot of corrections ~ g 2n to the instanton 
induced splitting of the ground state energy ~ e' 1 ^ 69 / g. The asymptotics of perturbative 
expansion of energy shift takes the form (see ||): 

Op- 1 / 6 ^ 2 3 n+2 r>' 
AE = — -=- (1 + ^ 9 2n El) , El = -(ln(6n) + 7 ) , (12) 

9VK n=l 71 
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where 7 = 0.5772... is the Eurlers constant. The perturbative series for the energy shift 
looks very similar with that for the ground state energy. The "best accuracy" one can 
reach after summing this series (which is not worse than the smallest term of the series) 
is negligibly small (~ exp(— l/2g 2 )). In order to reach the accuracy ~ exp(— l/3g 2 )/g 2 
in E + or 22_ one have to break the summation in (|12"D at the term for which g 2n E\ <C 
exp(— l/6g 2 ). The three - instanton contribution to AE, which generates the asymptotics 
(|j~2|) is considered in the Appendix (see (|38|)). 



3 Corrections to effective action 



Up to now we have taken into account only the leading term of interaction between instan- 
ton and antiinstanton AS* ~ e~ T /g 2 (fty. The further correction to the action of the pair 
will be considered in this section. 

The corrections of two types to the effective interaction are to be taken into account. 
The first is AS C ~ e~ 2T /g 2 - the correction to the classical instanton - antiinstanton 
interaction. The second correction ASg ~ e~ T appears after the Gaussian integration 
around the classical configuration. At first sight it seems that both these corrections do 
not contribute to the asymptotics of the perturbation theory. For example let us substitute 
the one loop correction ASg ~ e ~ T in the equation (0): 



E{3) = - / ^ exp { - 7 (1/3 - 2e ' T) } ( 1 + Be ") ■ 



(13) 



Here B ~ 1. Passing to the variable x = g 2 S we get 



E(g) 



1/3 -x 



1 + ^5(1/3 -x) 



(h (14) 



Tig 2 



The term proportional to B cancels in the singularity of the Borel function F(x) (|8|) and 
therefore does not contribute to the asymptotics of the perturbation theory. It is easy to 
show that correction to the classical action AS C ~ e~ 2T / g 2 also does not contribute to the 
asymptotics. 

Nevertheless as we will see below the corrections to the Gaussian integral and to the 
classical interaction appear which contain an additional large factor T: ASg ~ Te~ T and 
ASc ~ Te~ 2T / g 2 . Taking into account of these corrections leads to a new singularity of the 
Borel function SF ~ log(l/3 — x), although weaker than the pole (§). 

The quantum correction to the two instanton density is also to be taken into account: 

p = P/PA = ^exp(-25 / )(l+ Ag 2 ) . (15) 
ng 2 

Because of the correction Ag 2 does not depend on T, the value of A is twice the value of 
corresponding quantum correction to a single instanton density. The calculation of two 
loop correction even for single instanton is a complicated problem (see |TIJ), but one can 
avoid these calculations. The splitting of the ground state energy AE ([TJ) just equals to 
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the twice instanton density. The semiclassical method of finding the exponentially small 
splitting was given in |TTJ. The accuracy of semiclassical calculation allows to search for 
any of ~ g 2n corrections to small AE. The first term was found in ||: 

E\ = - = ~. (16) 
1 2 12 v ; 



3.1 The choice of classical trajectory 

In order to introduce the collective variable describing the relative positions of instanton 
and antiinstanton we multiply the functional integral by the Faddeev-Popov unit: 

1 = J 5Mx) - r f m f T ~ T) dT , (ifj\x) = / iP(t)x(t)dt . (17) 

Here T is the distance between pseudoparticles, r = r(T) - some function and if> — ip(t, T) 
- a vector in the functional space. The functional integration is to be performed over the 
hyperplanes orthogonal to the vector if). We would call the instanton - antiinstanton pair 
the configuration in the functional space which minimizes the action for a given value of 
collective variable r(T). The equation of motion now reads: 

^H = #, (18) 

OX 

where £ is a Lagrange multiplier, which allows to satisfy the condition (if)\xi-A) = t. 

In order to describe the instanton - antiinstanton pair one have to introduce the second 
collective variable to - the position of the center of the pair. But at least if ip{t) is symmetric 
under transformation t — to — > to — t, the second collective variable may be chosen in such 
a way, that it will not affect the classical configuration. 

There is a great freedom in the choice of if>(t,T) and r(T) (|I?|). The only restriction is 
that after we have separated the collective variables the functional integral over all the rest 
(quantum) variables must be well convergent. Technically the functional integration reduces 
to calculation of det S S ^Jr A ^ ~ the determinant of the second action variation in the field 
of the pair. Among the directions in the functional space there are two the "worse". These 
are the almost zero modes — the eigenf unctions of the operator 4-f having anomalously 
small eigenvalues. The vector if) from ( |i~7D should not be orthogonal to these modes. Under 
this conditions any function ip(t) would lead to a reliable instanton - antiinstanton type 
solution of the equation fll8P . In the following we suppose that ip(t) is a symmetric function, 
which differs from zero only near the centers of separate instanton and antiinstanton and 
goes to zero very fast both inside and outside the pair. 

3.2 Corrections to the classical trajectory and action 

Thus over a wide range between pseudoparticles if) ~ and eq. ( |18|) takes the form: 

- + x^ A = Sgxj_ A - 2g 2 x]_ A , (19) 
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Here 'dot' means the time derivative. Let us represent xi_a as a sum of instanton and 
antiinstanton field plus some correction: 

= xi{t + T/2) + x A {t - T/2) + 5x . (20) 

Here xi is an unperturbed instanton (f|) placed at t = —T/2 and xa is an antiinstanton 
placed at t = T/2. Solving eq. (|19|) iteratively in the intermediate region between instanton 
and antiinstanton, we get: 

e -T/2 -T 

xla = (e* + e - *) + {e 2t + e~ 2t + 6) + ... , 5x ~ Qe~ T / g . (21) 

9 9 

Outside the instanton - antiinstanton pair, at \t\ — T/2 ^> 1, 5x is neglegibly small. 
The action accurate to terms ~ Te~ 2T / g 2 reads: 

S[x T -a] ^ S[xj] + S[x A } + S[5x] + (22) 
+2([ dt 6S ^ Xl ^ (xa + Sx) + f dt^Y^-5x+( y x I + x A )Sx\ _ oo + x I x A \^ OD + 

\J — 00 OX J— 00 uX 

+ / dt(3g 2 x 2 x 2 A — 3gx A xi — QgxiXA^x) ) = 

= 1(1/3 - 2e- T - 12Te~ 2T + 0{e~ 2T )) . 

g 2 

Here we use (xj + x A )5x\o = and the equation of motion = 0. Corrections to 

(p2]) are of the order of 5S ~ e~ 2T / g 2 and depend on the choice of if) (|T7[) , (|T5D . But the 
contribution ~ Te~ 2T /g 2 in the formula fl22| ) is model independent. 



3.3 Corrections to determinant 

Now we have to find the determinant of the operator 

The correction AU equals to 

AU = -24e" T (24) 

over the wide range between instanton and antiinstanton (T/2 — \t\ ^> 1) and goes to zero 
outside the pair (|t| — T/2 ^> 1). Let us use the equality 

fet U) = exp ( tr '"( 1+ + i ))■ (25) 

where M = — + 1. Consider the logarithm in fl2"5D as the power series in the potential 
U. Because of the Green function of the operator M : G (t,t') = e"'*~* y2 falls down 
exponentially at large \t — t'\, the interference of U\ and t/2 i n (|25|) (as well as AU) may 
be neglected in the zeroth approximation. In this approximation the total determinant 
factorizes to the product of the single instanton determinants. There are two sources of 
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the corrections to (25) of the order of Te T . These are the contribution of the first order 
in AU and the interference in the lowest order of \J\ and £/ 2 : 



-fdtG(t,t)AU = -12Te~ T , 
dtdt'G 2 {t,t'){y I (t)U A (t') + Ui($)U A (t)) 



-12Te 



-T 



The full determinant up to model dependent corrections ~ 0(e T ) reads 



det 



W ) 



def 



-d 2 

a tt 



-4 + 1 



■\-2ATe~ 1 +0(e- 1 )) 



(26) 



(27) 



4 Results 



All the corrections (|T5|),(^),(^7J) should be combined in expression like (^) or ([?[). Now if 
one pass again to the variable x = g 2 S(T) = 1/3 — 2e~ T — 12Te~ 2T + ■ • ■ the Borel function 
F(x) (^J) may be found in a form: 

F(x) (y/^ + 3ln W* - x) + . . .) (l - ? j-g 2 ) . (28) 

The straightforward calculation leads now to the asymptotics (^). The first attempt to 
calculate numerically the 1/n correction to the asymptotics was made in [5]]. The authors of 
were the first who have calculated the leading (~ n\) asymptotics (0) of the perturbative 
expansion of ground state energy of the double well oscillator (|l|). However their numerical 
value for the 1/n correction was incorrect, being 103/36 instead of 53/18. In Table I. 
about 10 first terms of approximate expansion (^]) are compared with the exact values 0. 
Starting from n = 5 the error in E n never exceeds 12%. 

It is convenient to integrate by part the ~ ln(l/3 — x) term in F(x) (|28|). Now the 
expression for the "center of band" transforms to (compare with (|TTD): 



^ 2 f f i/3 e-'/g - e- 1 ^ dx e-Wfn(6) ^ 3"+^! 2n ) / 53 2 \ 

(29) 

The accuracy of this formula depends on N - the number of the terms of perturbative 
expansion, which were found exactly. In order to clarify this issue let us rewrite the integral 
in ( ]29T) through the principal value integral: 

dx " = P dx^- + e~W 2 (ln(3g 2 ) - j) . 

1/3 — x Jo 1/3 — x 

Now it is easy to subtract explicitly the N terms of the approximate expansion: 

p r°° e- x l g2 dx " 3 n+1 n!^ 2n _ p f°° {?>x) N e~ x l^ dx 
Jo 1/3 — XTcg 2 it Jo 1/3 — x irg 2 

Formula ( p9| ) reaches the best accuracy when one minimizes the value of fl3"0"|). This mini- 
mum occurs at \N — l/3g 2 | ~ 1. For such N the integral in r.h.s. of ([30]) (as well as the 
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minimum term of the perturbative expansion) is of the order of ~ 1/ ge~ l ^ a2 . The error 
in E c b ( |29"D is much smaller than the minimum term of the series 5E c b ~ g°e~ 1 ^ 3g . The 
integral in r.h.s. of (|30|) is extremely easy to estimate if N = l/3g 2 with some (large) value 
N. For such g the "best value" of E c b reads: 

Ea „ J> £ „ + 3 -^{,„ m + 7 + - ^ ) (31) 

Except for the new factor 1 — 53/18 N (or 1 — 53g 2 /6) our results ([31]), (|29|) agree with that 
presented in the literature (see ||). In Table II we compare the exact (numerical) values 
of E c b with that found from the approximate formula (|3T| ) for various values of N = l/3g 2 . 
It is seen that the error in (|3lD is of the order of ~ e"^ = e -1 / 39 . 



Asknowledgements Authors are grateful to V.L.Chernyak for helpful discussions. The 
work of P.G.S. was supported by International Science Foundation. 



5 Appendix 

In this Appendix we would like to consider the three instanton contribution to the asymp- 
totics of ~ g 2n corrections to the nonperturbative splitting AE = E- — E + of the Hamil- 
tonian ([!]) ground state energy. Previously this asymptotics was found within the usual 
quantum mechanical approach in || , and via analytical continuation of the multiinstanton 
contribution over g in ||. Although the trick with analytical continuation [[J leads to the 
correct asymptotics, the author of || have not shown, what the specific three instanton 
effects contribute to AE. 

Following H consider the functional integral with inverted (x — > 1/g — x) boundaries: 

e -E+L _ e -E_L = TrPe -HL = N f Dx(t)e- S ; L > 1 . (32) 

J x{ -L )+x{ L ) = 1 /g 

Only the n-instanton configurations with odd n contribute to TrPe~ HL . The single instan- 
ton contribution is trivial: 

(TrPe- HL ) =(2e- L/2 )L- . (33) 

V / JWg 



There are two distinct three-instanton contributions linear in L. First is the contribution 
from compact molecule - like three-instanton object having two interior degrees of freedom. 
The second contribution appears due to the effect of "excluded volume". Consider the pair 
of instanton and antiinstanton separated by interval T from each other and the single 
instanton outside the pair. While the integral over positions of the pair gives the factor L, 
the instanton can "run" only over the volume L — T. Thus the interference ~ LT should be 
considered in addition to the factorized ~ L 2 contribution of the pair and single instanton. 

Like it was done for the instanton pair (Q) we have to subtract the factorized three- 
instanton and instanton plus pair contributions from the contribution of three-instanton 
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molecule . After such regularization the molecule gives 

h = ( 2e -V2)_L_ Le -i/6 ff * J J dtdt'(e- s ^ - e~ s ^ - e' s ^ + e' 1 '^ , (34) 

where t,t' > - are the distances between the first and second and second and third 
pseudoparticles. S(t, if) = l/g 2 (l/3 — 2e~* — 2e~* + • ■ ■) is the three-instanton action (with 
1/Qg 2 excluded) and S p (t) = l/g 2 (l/3 — 2e~ t + - ■ ■) is the action of instanton - antiinstanton 
pair (|5|). Let us restrict the integration in (|3~4] ) by condition t,t' < T reg , where intermediate 
scale 1 <C T reg <C e 1//692 should cancel in the final result. Consider the integral 

K = f jdidi'(e- s ^-e^ 2 ) . (35) 



Here the integration should be performed over the range S(t,t') > 0. After passing to 
variables t and x = g 2 S(t,t') the integration over t may be done explicitly 

f l/3 e -x/g 2 _ e -l/3g 2 , ,1/3 _ x x 

K = l dx 1/3 2 K iZ r^) +r -)- (36) 

All the rest contributions in fl5^|) are proportional to T reg and naturally lead to cancellation 
of the proportional to T reg part of (j36|) . 



The effect of excluded volume leads to the linear in L three-instanton contribution of 
the form 

'> = (^J^-^- 1 "'" 1^ ^~T) (e-*CO _ e -W) . ,37) 

After the change of variables T — > x = g 2 S p (T) the integral in fl3~T| ) transforms to that of 
©• 

Finally the splitting of the ground state energy reads: 

(vrc/) 3 / 2 io 1/3 — a; v 2 ; v ; 



Expanding Zn(l/3 — x)/(l/3 — x) in powers of x one finds the series (12 ). 
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1.52 


1.07 


11 


4.6755 10 12 


1.44 


1.06 



Table I 



The exact values of the coefficients of perturbative expansion E n (see 0) and the ratio 
of approximate to exact value found in the leading (~ n\) approximation and with ~ 1/n 
correction taken into account. It is seen that starting from n = 5 the accuracy of the 
approximate formula (|2|) is not worse than 12% . 



N = l/(3<? 2 ) 


4 


6 


8 


10 


12 


rpexact 


0.4439 


0.43797 


0.44832 


0.459178 


0.467156 


rptheor 


0.4367 


0.44367 


0.44933 


0.459307 


0.467173 



Table II 



The exact numerical energy E^ act and the approximate value E^ eor found from the 
formula ([IT]) as a function of g 2 . As we expect the accuracy is ~ e~ N . 
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